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Abstract 

We show that all zero energy eigenstates of an arbitrary m-state quantum spin chain Hamil- 
tonian with nearest neighbor interaction in the bulk and single site boundary terms, which 
can also describe the dynamics of stochastic models, can be written as matrix product states. 
This means that the weights in these states can be expressed as expectation values in a Fock 
representation of an algebra generated by 2m operators fulfilling m 2 quadratic relations 
which are defined by the Hamiltonian. 
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1 Introduction 



A number of problems in many particle systems have been studied with the help of so called 
matrix product states. The idea of this technique is to express physical quantities such as 
ground state wave functions or correlation functions as products of operators acting on an 
auxiliary space and fulfilling algebraic relations defined by the Hamiltonian of the system. 
Introduced in the context of lattice animals [Q] the technique has been used to find ground 
states of various quantum spin chains [0-0. 

As it has been shown by Derrida, Evans, Hakim, and Pasquier || the steady state of the 
one species asymmetric exclusion process with open boundaries can be written as a matrix 
product state. Later works study this process in more detail |J-|§. Examples for asymmetric 
exclusion processes with two species were investigated in || [L0| with the help of the same 
technique. The algebra used in these studies is generated by as many operators as states a 
single site can take. It has the important property that it leads to recurrence relations for 
the steady state of systems of different lattice lengths. For other problems this algebra had 
to be generalized by enclosing additional operators. This was first done in the study of the 



dynamics of the asymmetric exclusion processe 0, [L2| . Another example is the algebra for 



the reaction-diffusion model studied in |]T3| which is generated by twice as many operators 
as the one from |J. In this Fock-algebra the recursion property is lost. The same is true for 
the algebra used in |14]]-|]17| for stochastic models with parallel updating. 



One of the questions coming up naturally in this context is the following: To which kind 
of problems can the matrix technique be applied? In other words: Does a Hamiltonian, 
which describes either a quantum spin system or a stochastic process, need to have any 
particular property in order to have matrix product eigenstates? In this paper we will 
prove the following proposition: Any zero energy eigenstate of a Hamiltonian with nearest 
neighbor interaction in the bulk and single site boundary terms can be written as a matrix 
product state with respect to the Fock-algebra which is that of Ref. /|73|/. The 

technique of matrix product states was called matrix product ansatz in the literature. What 
we will see is that this technique, in the form of ( ^.2| )-( P73| ), is not an ansatz but rather an 
identical reformulation of the eigenvector equation for the zero energy eigenstate. 

The bulk of the paper is organized as follows: In Section 2 we define the class of Hamilto- 
nians and in Section 3 we give the proof of the proposition. We conclude with some remarks 
on matrix product techniques in Section 4. 



2 Definition of a class of models 

The Hamiltonian we are going to consider in this paper is of the form 

JV-l 

H = h^ + h iJ+i + h(R) (2-1) 
i=i 

with 



ftW = fc(0 (g, , frW = ® (2.2) 



1 



and h 



(2.3) 



where / is the m x m-identity matrix, /i is an m 2 x m 2 - matrix describing a two site interactions 
in the bulk, and , are m x m-matrices defining single-site boundary terms. This 
kind of Hamiltonian appears as well in the study of one-dimensional stochastic systems as 
in the study of quantum spin chains. We will now describe both applications in more detail. 

Stochastic model: Let us consider a one-dimensional lattice with N sites, each of which 
can be in either of m states. A configuration on the lattice is completely defined by the set 
of occupation numbers {s^} = s±, S2, ■ ■ ■ , sjv with Si = 1, m V i = 1, N . The system 
evolves stochastically. During an infinitesimal time step dt its configuration {s^ can change 
to a configuration {s^} with a probability r({sj}, {s'^dt where the r({si}, {s-}) are referred 
to as rates. This process can be described in terms of a rate equation which reads 



(2.4) 



where P({si},t) = P{ si, S2, ■ ■ ■ , sjvi t) is the probability of finding the configuration si, S2, ■ ■ ■ , sjy 
at time t . Throughout this paper we restrict ourselves to dynamics where the configuration 
can change only at two adjacent sites at a time and the rate for such a change depends only 
on these two sites. The rates are assumed to be independent of the position in the bulk of 
system. At the boundaries, i.e. at sites 1 and N, we assume additional processes to take 
place. 



(2.5) 



It is convenient to introduce a vector notation ]T3| by writing 
\P(t)) = J2 P(sus 2 , 

{Si} 



.,S N ,t) I Si ) ® |s 2 ) ® ... ® \s N ) 



with 



\k) 



(2.6) 



\0/ <- m 

In terms of these vectors the rate equation reads 



dt\P(t) 



-H\P(t) 



(2.7) 



where H is an m x m -matrix which is defined in terms of the rates. We call H the 



Hamiltonian here. For the processes we are studying it has the structure of ( p.l| )-( pT3|) . For 
stochastic models the matrices H have a particular property: They have to have vanishing 
column sums because the total probability has to be conserved. This implies that there is a 



2 



zero energy eigenstate for every lattice length N which is just the steady state of the system. 
It is denoted by \Pn ) and obeys the relation 

H\P N ) =0 . (2.8) 

In Section 3 we will discuss a representation of these states. 



Quantum spin chain: Consider a chain of length N with a spin s-particle sitting on 
each site. Suppose there is an interaction between adjacent particles and some surface fields 
acting on sites 1 and N . Then the system dynamics is described by a Schrodinger equation 
with a Hamiltonian of type ( p.l| )-( pT3|) if we chose a spin s-representation and m = 2s + 1 . 
The Hamiltonian generally does not have vanishing columns sums like in the stochastic case, 
though it has to be hermitian. This condition does not effect the construction of matrix 
product states below. Although the physical meaning of the Schrodinger equation is different 
from the meaning of the rate equation ( |2.7|) it has the same mathematical structure. And 
just as for the stochastic model we are often interested in the ground state of the Hamiltonian 
H. 



3 Matrix product states 



Let us now turn to the matrix product states 0-|13|. We introduce an auxiliary vector 
space V a and we define 2m operators D s and X s with s = 1, 2, m acting on V a as well as 
two vectors \W > and < V\ in V a . ( Vectors in the auxiliary space are denoted by |... > 
in contrast to the vectors in the configuration space which are denoted by |... ) .) Next we 
define a vector \P N ) as 

f D 1 \® N 
D 2 



\Pn) =< W\ 



w > 



(3.1) 



\D m ) 



where ® stands for the direct product in the configuration space, so that the above equation 
reads in terms of its components: P^{s\, s 2 , Sjy) =< W\D S1 D S2 ...D SN _ 1 D SN \V > . We 
call states of type \Pn ) matrix product states. 



Let H be a Hamiltonian of type (|2.1| )- (|2.3j ) which has a zero energy eigenstate for all 
lattice lengths iV (a Hamiltonian describing a stochastic model always has this property). 
Then we can state the following proposition: 



(i) If D s , X s and \W > , < V\ fulfill the following relations 



h 



( D x \ 

D 2 



D 2 



x 2 



( D x \ 
D 2 



D 2 



X 2 



(3.2) 



\D m j 



\D m j 



\X m J 



\D m ) 



\D m j 



\X m J 
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< W\ h® 



(D 1 \ 
D 2 



\D m J 



< W\ 



x 2 



\X m J 



and h^ r ' 



f D 1 \ 

D 2 



\D m J 



\V >- 



(x 1 \ 

x 2 



\X m J 



w > 



(3.3) 



then the vector \Pn ) defined by Eq. solves the equation for the zero energy eigenstate, 

i.e., 

H\P N ) =0 . (3.4) 

(ii) For any vector |Pjv ) solving H |Pjv ) = one can find operators D s , X s (s = 1, 2, m) 
and vectors < W\ , \V > in some space V a , such that \Pn ) can be represented as a matrix 
product state \Pn ) defined by Eq. (ffi 1\) and relations (\3.3j) and (\3.c\) are fulfilled in V a . 



Before we come to the proof of the statement we stress again that, in contrast to the 
algebra defined in the relations (|3.2|) - (|3.3|) do not lead to recurrence relations for 

expectation values of products of different length in the operators D s and X s . 

The proof of (i) is based on a site by site cancellation of terms when H is applied on 
| Pjv ) • The mechanism is basically the same as the one worked out in Ref . f5[ . In Appendix 
A we explain it in detail. 

In order to prove (ii) we construct a representation of operators and vectors fulfilling the 
Fock-algebra ( |3.2|) -( [3~3D using the eigenvectors \Pn ) of systems with lengths N = 1,2, ... . 
Let us define for any number M = 1,2,... an m M -dimensional space Vm with a set of 
orthogonal basis vectors \si, s 2 , sm > (^i = l,2,...,m) as well as the one-dimensional 
space Vq with the basis vector | > . We define the space V a as the direct sum of all Vm with 
M = 0,1,2,... . The above basis vectors may be written as infinite column vectors with a 
1 at position 1 + s\ + ms 2 + m 2 ss + .... + m M ~ 1 SM and a at all other positions, and their 
transposed < si,S2,...| can be written as the corresponding row vectors. Next we define 
operators D s and X s as well as vectors < W\ and \V > by means of their action on all the 
|si, s 2 , ... > which provides a matrix representation if we write the \si,s 2 , ... > as columns. 
The D s we define by 

D si \s 2 ,s 3 , ...,s N > = \s 1 ,s 2 ,s 3 , ...,s N > (3.5) 

for N=l,2, ... and s, = l,2,...,m . 

The X s act on the basis vectors of V/v_i , i.e. on the \si,s 2 , sat_i > , as 

/ N— 1 numbers \ 





IM, 


...,i> 


fxA 




...,2> 


x 2 
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\X m J 







\m, m, m — 1 > 
\ \m, m, ...,m > J 
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[hi,2 + h 2<3 + ... + h N _ 1<N + h (R ^} 



( D x \ 
D 2 



\D m ) 



> 



(3.6) 



for N = 1,2, 



where /i^i+i and ft,^ act on the N-fold tensor product in Eq. ( |3.6D according to their 
definitions (2.2) and ( j2.3| ). Note that each choice of iV gives a different set of basis vectors 
in the column on the l.h.s., and taking successively N — 1, 2, ... all basis vectors will occur in 
this column exactly once. Hence, the above equation defines the action of each X s on each 
basis vector in a consistent way. The vectors \V > and < W\ we define by 



\V > 

< W\s 1 ,s 2 , > 



| > 

Pn(si,s 2 , .-.,s n ) 

for N = 1,2,... and s { = 1,2, 



(3.7) 
(3.8) 



, m 



where Pjv(si, s 2 , Sjv) are the components of |Pjv ) • The above definitions fix the D s , X s 
and the < W\ , \V > completely up to a constant < W\ > . This constant can be fixed 
arbitrarily, since it does not enter into any result. 

We have to prove now that these operators and vectors lead back to |Pjy ) by means of 
Eq. (|3.1|) , i.e., that |Pjy ) is equal to |P/y ) , and that they obey the Fock-algebra (|3.2|) - (|3.3|) . 
Let us begin the first proof by rewriting ( |3.5| ) as 



D Sl D S2 ...D SN \ >= \s u s 2 , s N > 
for iV=l,2, ... and Sj = l,2,...,m . 



(3.9) 



(The above relation gives a more intuitive meaning of the representation we have chosen: 
A basis vector \si, s 2 , sm > is created by applying the sequence D Sl D S2 ...D SM on | > .) 
Furthermore, using (|3.8|) we get: 



< W\D S1 D S2 ...D SN \ >= P N (s 1 ,s 2 , ....,sjv) 



(3.10) 



for N = 1,2,... and st = 1,2, 



m 



Rewriting (|3.10|) as a direct product of vectors in the configuration space and using | >= 
\V > yields 



fDi\ 
D 2 



< w\ 



w > 



IP 



N 



(3.11) 



\D m j 

which is what we were to prove 



for AT = 1,2, 
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In order to show that above representation fulfills the algebra ( |3.2|) -( |3~3"1) we rewrite the 

D 2 

definition ( |3.6|) using Eq. ( |3.9|) and multiply the column vector 



from the left: 



\D m ) 



D 2 



(Xx\ 

x 2 



(Di\ 
D 2 



®(JV-1— n) 



> 



\D m ) \X m ) \D m j 



[h n +l,n+2 + ^n+2,n+3 + ••• + ^JV-1,JV + h^} 



D 2 



> 



(3.12) 



\D m ) 

for n = 0,1,..., N- 1 and AT =1,2,... . 



The right equation of (^3) is nothing but the case n — , iV = 1 of above relation. 
Consequently it holds. Next we use Eq. ( p.l2| ) to compute 



Xo 



(D 1 \ 
Do 



®(iV-l) 



(Di\ 
D 2 



> 



\X m ) \D m J 

fD 1 \ 9N 
D 2 



Xo 



( D x \ 

Do 



®(JV-2) 



> 



\D m y \x m / \D m y 



> 



(3.13) 



i.e., 



f (Xi\ 

Xo 



\D m ) 



Do 



(D 1 \ 

Do 



Xo 



h 



\X m J \D m J \D m J \X m J 

(D,\ m ~ 2) 
D 2 

| > = 

\D m ) 
for iV = 2,3,... . 



Do 



(D 1 \ 

Do 



\D m j \D m j J 



(3.14) 
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Since the products of the D s span the whole space V a , Eq. (|3.2p must hold in V a . What 
remains to be proven is the left equation of ( |3.3| ). We multiply Eq. ( |3.12| ) for n = by < W\ 
from the left and use Eq. ( |3.11| ): 



< W\ 



\X m J 



(D x \ 
Do 



\D m ) 



i(N-l) 



> 



[h,2 + h 2 , 3 + ... + h N . ltN + \P N ) . 



(3.15) 



Writing now the Hamiltonian of an iV-site system as H = h^ L ' + h lt2 + h 2t3 + ... + /iAr_i,jv + /i^ R - ) 
and using H\P^ ) = , the above relation can be written as: 



< w\ 



® 



\X m J 



\D m j 



i(N-l) 



> 



-ftW IP, 



A' 



(3.16) 



Using again Eq. ( |3.11| ) yields 







(DA 




(XA ) 


(DA 






D 2 




x 2 


D 2 


<w\ < 






+ 


. 






- 


\D m ) 




\X m J 


\D m ) 



®(JV-1) 



> = 



(3.17) 



for N = 1,2, 



Consequently the left equation of ( |3.3|) holds in V a . This completes the proof that the 
matrices defined by Eq. fl3.5|)-(|3~8D represent the algebra (|3.2|)-(|3~3|) and hence the proof of 
proposition (ii). 



Let us add some remarks on the proof. We proved that the Fock algebra (|3.2| )- (|3.3| ) has 



a nontrivial representation if an eigenstate with zero energy exists for all lattice lengths. 
The construction of the matrices obeying ( |3.2| ) can be done without any demand for special 
properties of H. (An associative algebra generated by 2m generators with m 2 quadratic 
relations is always a well defined mathematical object.) The special property of H, i.e. the 
existence of zero energy eigenstates \Pn{s\, s 2 , sjv)), was only needed for the construction 
of a non trivial scalar product involving the vectors < W\ and \V >. More precisely, the 
vectors \Pn(si, s 2 , Sjv)) only enter the definition of the vector < W\ and the equation 
which determines < W\ is just the equation for the zero energy eigenstate. 



Due to the lack of recurrence relations the scalar product in fl3.2| )- fl3~3"| ) is not completely 
fixed and can be chosen independently in each vector space Vm- For that reason it would be 
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sufficient to require that at least for one lattice length TV the Hamiltomian has a zero energy 
eigenstate. In that case < W\ is nontrivial only on the subspace Vjy. 

One may define a finite dimensional V a as the direct sum over all Vm with M < N max for 
some number N max . A matrix representation of the algebra ( p.2| )-( pT3|) can then be found 
be assuming Eq. (|3^)-(^8|) for N = 1,2, N max as well as D s \si, s 2 , s Nmax >= and 
X s \si, S2, sjsr max >= . Then Eq. fl3.1|) provides a matrix product representation for all 
eigenvectors \Pn ) with N < N max . 

It is of course also possible to express eigenstate of non zero eigenvalues as matrix product 
state by adding an appropriate shift to the Hamiltonian. This is done by replacing h, 
and hS r > in Eqs. ( |3.2| )-( fr3|) by h — e, — and — e^ r \ respectively. In this case the 
Hamiltonian H in Eq. (g§ has to be replaced by H-E(N) with E(N) = S + e (r) + (N -l)e 
if N is the lattice length. The vector \Pn) is then an eigenstate with energy E(N). However, 
if e W,e (r) and e are independent of N, Eq. ( J3.4| ) will have a non trivial solution only for 
a finite number of lattice lengths N and < W\ will therefore be non trivial only on the 
corresponding subspaces Vjv of V a . 



4 Concluding Remarks 

We proved that any zero energy eigenstate of a Hamiltonian of the form can be written 
as a matrix product state with respect to the Fock representation of the algebra (|3.2| ) -(|3.3| ) . 
The proof has been done by showing that the operators as well as the scalar product can 
be defined in a non trivial way. The equations defining the scalar product are exactly the 
equations for the zero energy eigenstate. This shows that by application of the algebra on 
an abstract level, i.e. by using only the relations (p.2|)-( p73|) , one can not gain any insight 
into the form of the eigenstates. One ends up with nothing but a reformulation of the 
equation for the zero energy eigenstate. Therefore the so called matrix product ansatz, i.e. 
the application of the Fock-algebra (|3.2|) - (|3.3| ) to an eigenvector problem, is not really an 
ansatz, it is an identity. The eigenvector problem is not transformed into another problem, it 
remains unchanged. However, the situation becomes different if one considers a more special 
form of the algebra as done in ]5|-|JlU|j where the operators X s are replaced by numbers. 

We checked that our proposition is also true for the Fock-algebra describing stochastic 



models with parallel updating []14[ . 
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A Appendix: Proof of eigenvalue equation for matrix 
product state 



Following the line of Ref. || we prove proposition (i), i.e., we show that the states ( |3.1| ) 
solve the eigenvalue equation (|3.4j ) if (|3.2|) -(|3~3l) are fulfilled. Using the definition (|2TT| ) we 
write 



(N-l) 

i=i 

+ I®*"" 1 ) (g) fcM 



Next we define a state |Pjy(fc) ) as 



|Pjv(fc) ) = < W\ 



D 2 



\D m ) 



i(k-i) 



fx x \ 

x 2 



\X m J 



( D x \ 
D 2 



i(N-k) 



w > 



\D m ) 



Applying now the terms occurring in Eq. ([A.1Q on the state 



\Pn) =<W\ 



D 2 



\D m ) 



w > 



yields 



h (i) \p N j 

J®*"- 1 ) ® fcM \P N ) 



-_\Pn(1)) 

\P_nU)) ~ \PnU + 1)) 
\Pn(N) ) 



(A.l) 



(A.2) 



(A.3) 



(A.4) 
(A.5) 
(A.6) 



where the algebra ( |3T2|) - (|3.3| ) was used. Applying now the full operator H (see Eq. ( |A.1| )) 

on the state |P/v ) we get 



N 



H IP, 



N 



-IPjv(I)) + J2(\ P n(j)) ~ \P N (j + l) 




\Pn{N) 



(A.7) 



which is what we were to show. 
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